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Abstract 

We prove that for each k > 0, the probability that a root vertex in a random planar graph 
has degree k tends to a computable constant dk, so that the expected number of vertices 
of degree k is asymptotically d^n, and moreover that dk — 1. The proof uses the tools 
developed by Gimenez and Noy in their solution to the problem of the asymptotic enumeration 
of planar graphs, and is based on a detailed analysis of the generating functions involved in 
counting planar graphs. However, in order to keep track of the degree of the root, new 
technical difficulties arise. We obtain explicit, although quite involved expressions, for the 
coefficients in the singular expansions of the generating functions of interest, which allow us 
to use transfer theorems in order to get an explicit expression for the probability generating 
function p(w) = ^ fe dkW k . From this we can compute the dk to any degree of accuracy, and 
derive the asymptotic estimate dk ~ c • k~ 1 ^ 2 q k for large values of k, where q » 0.67 is a 
constant defined analytically. 

1 Introduction 

In this paper all graphs are simple and labelled with labels {1,2, . . . ,n}. As usual, a graph is 
planar if it can be embedded in the plane without edge crossings. A planar graph together with 
a particular embedding in the plane is called a map. There is a rich theory of counting maps, 
and part of it is needed later. However, in this paper we consider planar graphs as combinatorial 
objects, regardless of how many non-equivalent topological cmbcddings they may have. 

Random planar graphs were introduced by Denisc, Wasconcellos and Welsh [10], and since 
then they have been widely studied. Let us first recall the probability model. Let Q n be the family 
of (labelled) planar graphs with n vertices. A random planar graph lZ n is a graph drawn from Q n 
with the uniform distribution, that is, all planar graphs with n vertices have the same probability 
of being chosen. As opposed to the classical Erdos-Renyi model, we cannot produce a random 
planar graph by drawing edges independently. In fact, our analysis of random planar graphs relies 
on exact and asymptotic counting. 

Several natural parameters defined on lZ n have been studied, starting with the number of 
edges, which is probably the most basic one. Partial results where obtained in [TDl [17l [28] E], 
until it was shown by Gimenez and Noy [19] that the number of edges in random planar graphs 
obeys asymptotically a normal limit law with linear expectation and variance. The expectation is 
asymptotically nn, where k ~ 2.21326 is a well-defined analytic constant. This implies that the 
average degree of the vertices is 2k f=a 4.42652. 

McDiarmid, Steger and Welsh [25] showed that with high probability a planar graph has a 
linear number of vertices of degree k, for each k > 1. Our main result is that for each k > 1, the 
expected number of vertices of degree k is asymptotically dkn, for computable constants dk > 0. 
This is equivalent to saying that the probability that a fixed vertex, say vertex 1, has degree k 
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tends to a limit dk as n goes to infinity. In Theorem 16.131 we show that this limit exists and we 
give an explicit expression for the probability generating function 



from which the coefficients dk can be computed to any degree of accuracy. Moreover, we show 
that p(w) is indeed a probability generating function, that is, ^ dk = 1. 

The proof is based on a detailed analysis of the generating functions involved in counting planar 
graphs, as developed in [19] . where the long standing problem of estimating the number of planar 
graphs was solved. However, in this case we need to keep track of the degree of a root vertex, and 
this makes the analysis considerably more difficult. 

Here is a sketch of the paper. We start with some preliminaries, including the fact that for the 
degree distribution it is enough to consider connected planar graphs, and that d = 0. Then we 
obtain the degree distribution in simpler families of planar graphs: outcrplanar graphs (Section [3]) 
and series-parallel graphs (Scction[4|. We recall that a graph is series-parallel if it does not contain 
the complete graph K± as a minor; equivalently, if it does not contain a subdivision of Since 
both K§ and ^3,3 contain a subdivision of K4, by Kuratowski's theorem a series-parallel graph 
is planar. An outcrplanar graph is a planar graph that can be embedded in the plane so that all 
vertices are incident to the outer face. They are characterized as those graphs not containing a 
minor isomorphic to (or a subdivision of) either K4 or ^2,3- These results are interesting on their 
own and pave the way to the more complex analysis of general planar graphs. We remark that the 
degree distribution in these simpler cases has been obtained independently in [5J [5] using different 
techniques. 

In Section Owe compute the generating function of 3-connected maps taking into account the 
degree of the root, which is an essential piece in proving the main result. We rely on a classical 
bijection between rooted maps and rooted quadrangulations [9l |26|, and again the main difficulty 
is to keep track of the root degree. 

The task is completed in Section[6l which contains the analysis for planar graphs. First we have 
to obtain a closed form for the generating function B'(x, y, w) of rooted 2-connectcd planar graphs, 
where x marks vertices, y edges, and w the degree of the root: the main problem we encounter 
here is solving a differential equation involving algebraic functions and other functions defined 
implicitly. The second step is to obtain singular expansions of the various generating functions 
near their dominant singularities; this is particularly demanding, as the coefficients of the singular 
expansions are rather complex expressions. Finally, using a technical lemma on singularity analysis 
and composition of singular expansions, we are able to work out the asymptotics for the generating 
function C*(x,y,w) of rooted connected planar graphs, and from this the probability generating 
function can be computed exactly. We also compute the degree distribution for 3-connected and 
2-connected planar graphs. Finally in Section [7] we show that there exists a computable degree 
distribution for planar graphs with a given edge density or, equivalently, given average degree. 

For each of the three families studied we obtain an explicit expression, of increasing complexity, 
for the probability generating function p(w) = J2k>i dkW k . Theorems 13.31 B~5l and 16. 131 give the 
exact expressions in each case. We remark that the expression we obtain for p(w) in the planar 
case is quite involved and needs several pages to write it down. However, the functions involved 
are elementary and computations can be performed with the help of Maple. 

The following table shows the approximate values of the probabilities dk for small values of k, 
which are obtained by extracting coefficients in the power series p(w), and can be computed to 
any degree of accuracy. 

We also determine the asymptotic behaviour for large k, and the result we obtain in each case 
is a geometric distribution modified by a suitable subexponential term. We perform the analysis 
for connected and 2-connected graphs, and also for 3-connected graphs in the planar case. Table 
[2] contains a summary of the main results from sections [3j 0] and [6l It is worth noticing that 
the shape of the asymptotic estimates for planar graphs agrees with the general pattern for the 
degree distribution in several classes of maps, where maps are counted according to the number 
of edges [22] . 
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Table 1: Degree distribution for small degrees. 
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Table 2: Asymptotic estimates of dk for large k. The constants c resp. c! and q in each case are 
defined analytically. The two approximate values in the last row are exactly the same constant. 

As a final remark, let us mention that in a companion paper [12] , we prove a central limit 
theorem for the number of vertices of degree k in outerplanar and series-parallel graphs, together 
with strong concentration results. It remains an open problem to show that this is also the case 
for planar graphs. Our results in the present paper show that the degree distribution exists and 
moreover can be computed explicitly. 



2 Preliminaries 

For background on generating functions associated to planar graphs, we refer to [H] and [7J, 
and to [57] for a less technical description. For background on singularity analysis of generating 
functions, we refer to the Q2] and to the forthcoming book by Flajolet and Sedgewick [H] , 

For each class of graphs under consideration, c n and b n denote, respectively, the number of 
connected and 2-connected graphs on n vertices. For the three graphs classes under consideration, 
outcrplanar, series-parallel, and planar, we have both for c n and b n estimates of the form 

c-n- a p- n n\, (1) 

where c, a and p are suitable constants pj 119) . For outcrplanar and series-parallel graphs we have 
a = —5/2, whereas for planar graphs a = —7/2. A general methodology for graph enumeration 
explaining these critical exponents has been developed in [20] . 

We introduce the exponential generating functions C(x) = c n x n /nl and B(x) = b n x n /n\. 
Let Cfc be the exponential generating function (GF for short) for rooted connected graphs, where 
the root bears no label and has degree k; that is, the coefficient [x n /n\]Ck{x) equals the number 
of rooted connected graphs with n + 1 vertices, in which the root has no label and has degree k. 
Analogously we define Bk for 2-connected graphs. Also, let 

B%x,w)=J2 B k{x)w k , C'{x,w) = Y. C ^ x ) wk - 

k>2 k>0 

A basic property shared by the classes of outerplanar, series-parallel and planar graphs is that 
a connected graph G is in the class if and only if the 2-connected components of G are also in the 
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class. As shown in |19j . this implies the basic equation 

C'{x) = e B '^ c "^ 

between univariate GFs. If we introduce the degree of the root, then the equation becomes 

C'(x,w) = e B ' (xC ' {x) < w) . (2) 

The reason is that only the 2-connected components containing the root vertex contribute to its 
degree. 

Our goal in each case is to estimate [x n \Ck(x), since the limit probability that a given fixed 
vertex has degree k is equal to 

[x n ]C k (x) 

4 = »tTO' (3) 

Notice that in the denominator we have the coefficient of C'(x), corresponding to vertex rooted 
graphs in which the root bears no label, in agreement with the definition of C k {x). 

A first observation is that the asymptotic degree distribution is the same for connected members 
of a class than for all members in the class. Let G(x) be the GF for all members in the class, and 
let Gk(x) be the GF of all rooted graphs in the class, where the root has degree k. Then we have 

G(x)=e ci - X \ G k (x)^C k (x)e c ^. 

The first equation is standard, and in the second equation the factor C k {x) corresponds to the 
connected component containing the root, and the second factor to the remaining components. 
The functions G{x) and C(x) have the same dominant singularity. Given the singular expansions 
of G(x) and C(x) at the dominant singularity in each of the cases under consideration, it follows 
that 

[x n ]G k (x) = Um [x n }C k (x) 
[x n ]G'(x) [x n ]C'(x) ' 

Hence, in each case we only need to determine the degree distribution for connected graphs. A 
more intuitive explanation is that the largest component in random planar graphs eats up almost 
everything: the expected number of vertices not in the largest component is constant |23j . 

Another observation is that do = and d\ = p, where p is the constant appearing in the 
estimate ^ for c n ; as we are going to sec, p is the radius of convergence of C{x). Indeed, there 
are no vertices of degree zero in a connected graph, and the number of connected graphs in which 
the root has degree one is n(n — l)c n _i / 'nc n ~ p. 

The general approach we use for computing the d k is the following. Let f{x) = xC'(x) and let 
H(z) = e B '( z ' w \ where w is considered as a parameter. Let also p be the radius of convergence 
of C(x), which is the same as that of f(x). According to ([2]) we have to estimate [x n ]H(f(x)), 
and this will depend on the behaviour of H(z) at z = f(p). In the outerplanar and series-parallel 
cases, H{z) turns out be analytic at f(p), whereas in the planar case we have a critical composition 
scheme, that is, the dominant singularity of H(z) is precisely f(p). This is a fundamental difference 
and we have to use different tools accordingly. Another difference is that B* is much more difficult 
to determine for planar graphs. 

Finally we comment on an asymptotic method that we apply several times. Suppose that 
f( z ) = 12n>o a ™ z ™ i s the power series representation of an analytic function and p > is the 
radius of convergence of f(z). We say that f{z) is analytic in a A-region if f(z) can be analytically 
continued to a region of the form 

A = {zG<C: \z\ <p + V, \axg(z-p)\>6}, (4) 

for some r\ > and < 6 < J. 

If we know that |/(z)| < C ■ |1 - z/p\- a for z G A, then it follows that \a n \ < C ■ p^n 01 ' 1 
for some C > that depends on C,a,rj, and 9; see [13] . In particular, if we know that f(z) is 
analytic in a A-region and has a local representation of the form 

f(z) = A + A 2 Z 2 + A 3 Z 3 + O(Z 4 ) 1 (5) 
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where Z = sj\ - z/p, then it follows that \f(z) - A - A 2 Z 2 - A 3 Z 3 \ < C- |1 - z/p\ 2 for zeA. 
As a consequence 

a n = ^±p-n-^ + 0(p- n n- 3 ). 

In fact, we focus mainly on the derivation of local expansions of the form j5]). The analytic 
continuation to a A-rcgion is usually easy to establish. We either have explicit equations in known 
functions or implicit equations where we can continue analytically with the help of the implicit 
function theorem. 

As a key example, we consider a function y = y{z) that has a power series representation 
at Zq = and that satisfies an analytic functional equation $(y, z) = 0. Suppose that we have 
y(zo) = yo (so that &{yo, z ) = 0) and § y (yo, z ) ¥^ 0- Then the implicit function theorem implies 
that y{z) can be extended analytically to a neighbourhood of z = z$. In particular it follows 
that y(z) cannot be singular at z = z$. On the other hand if we know that there exists zq and 
2/o = y(zo) with 

$(y ,zo)=0 and <Z> y (y , z ) = 

and the conditions 

®z(yo, z ) ^ and ® yy (yo, z ) ^ 0, 
then Zq is a singularity of y(z) and there is a local expansion of the form 

y{z) = Y + Y 1 Z + Y 1 Z 2 + ■■■ , 

where Z = yjl- z/z , Y = y and Y\ = -^/2z $ z (y , z )/^ yy (y 0l z Q ); see [TT]- In our appli- 
cations it is usually easy to show that <& y (y(z), z) =/= for \z\ < zq and z ^= Zq. Hence, in this 
situation z = zq is the only singularity on the boundary of the circle \z\ < zq and y(z) can be 
analytically continued to a A-region. 



3 Outerplanar graphs 

In this section C(x) and B(x) now denote the GFs of connected and 2-connected, respectively, 
outerplanar graphs. We start by recalling some results from [7]. From the equivalence between 
rooted 2-connected outerplanar graphs and polygon dissections where the vertices are labelled 
1,2, ... ,n in clockwise order (see Section 5 in [7J for details), we have the explicit expression 



_.. , l + 5x— -v/l — 62; + ; 
B ( x > = S 



The radius of convergence of B{x) is 3 — 2\/2, the smallest positive root of 1 — 6x + x 2 = 0. The 
radius of convergence of C(x) is p = ^p(r), where ip(u) = ue~ B W, and r is the unique positive 
root of ip'(u) — 0. Notice that t satisfies tB"(t) = 1. The approximate values are r w 0.17076 and 
p f=a 0.13659. We also need the fact that ip is the functional inverse of xC'(x), so that r = pC'(p). 
Let 

1 + x - \/l - 6x + x 2 
D(x) = (6) 

and let D^{x) = x(2D(x) — x)^ 1 (D/. is the ordinary GF for polygon dissections in which the 
root vertex has degree k). Then we have 

B k = ^D k , k>2, B 1= x. 
By summing a geometric series we have an explicit expression for B* , namely 



00 



Our goal is to analyze B'(x,w) and C'(x,w) = exp(B* (xC (x) , w)) . For this we need the 
following technical lemma. 



5 



Lemma 3.1. Let f(x) = ^™>n o, n x n /n\ denote the exponential generating function of a sequence 
a n of non-negative real numbers and assume that f{x) has exactly one dominating square-root 
singularity at x = p of the form 



f(x) = g{x) - h(x)y/l- x/p, 

where g(x) and h{x) are analytic at x = p and f(x) has an analytic continuation to the region 
{x £ C : \x\ <ple}\{i£l:i>p} for some e > 0. Further, let H(x, z) denote a function that 
is analytic for \x\ < p + e and \z\ < f(p) + s such that H z (p, f(p)) 7^ 0. Then the function 

f H (x) = H(x,f(x)) 

has a power series expansion fn(x) = J2 n >a °nX n /n\ and the coefficients b n satisfy 

lun^=H z (p,f(p)), (8) 

n— >oo a„ 



Proof. From f(x) = g(x) — h(x)^/l — x/p it follows from singularity analysis |14| that the sequence 
a n is given asymptotically by 

«n Hp) „ 3/2 

n\ ~ 20F 9 ■ 
Since H is analytic at f(p), it has a Taylor series 

H(x, z) = H(p, f{p)) + H z (p, f(p))(z - f(p)) + H x (p, f( P ))(x -p) + ... 

The function fn(x) has also a square-root singularity at x = p with a singular expansion, obtained 
by composing the analytic expansion of H (x, z) with the singular expansion of f(x), namely 



f H (x) = H(p, f(p)) - H z (p, f(p))h(p)Jl -- p + 0(\1 - x/p\). 

Consequently, the coefficients b n can be estimated as 

bn H z (pJ{p))h(p) _ n _ 3/2 
n\ ~ 2^ 9 n ' 

and © follows. □ 

We are ready for obtaining the degree distribution of two-connected outerplanar graphs and 
connected outerplanar graphs. Both results have been obtained independently in [3 [6], and our 
respective results agree. 

Theorem 3.2. Let dk be the limit probability that a vertex of a two-connected outerplanar graph 
has degree k. then 



*»> - g - ir^ra - g 2(3 - 2 ^< 4 - - '» 

Moreover p(l) = I, so that the dk are indeed a probability distribution. 

Proof. Since B'(x 1 1) = B'(x) and D(x) = 2B'(x) — x we can represent B*(x, w) as 

xw 2 AB'(x)-?>x 
B'(x,w) = xw -\ ; — ^— • 

v ' ; 2 l-(4B'(x) -3x)w 

Hence, by applying Lemma |3 . 1 1 with f(x) = B'{x) and 

xw 2 Az — 3x 

H(x, z) = xw + 



k w k . 



2 l-(4z- 3x)i 
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we obtain 

,,2 



p(w) = lim 



[x n ]B*{x,w) 2(3-2^2)^ 



n^oo [x n ]B'(x) (1 - (y/2 - l)w) 2 ' 

Note that p = 3 — 2\/2 and that w is considered here as an additional (complex) parameter. □ 

Theorem 3.3. Let dk be the limit probability that a vertex of a connected outerplanar graph has 
degree k. then 



k>l 



x=pC'(p) i 



where B* is given by Equations (0|) and 

Moreover p{\) — 1, so that the dk are indeed a probability distribution and we have asymptoti- 
cally, as k — > oo 

where c x w 0.667187, c 2 « 0.947130, and g = 2D(r) - r « 0.3808138. 
Proof. We have 

^ C fe (a> fc = C*(a;, «;) = e B '^ c '^' w \ 
k 

The radius of convergence 3 - 2^2 of B(x) is larger than pC'(p) = r « 0.17076. Hence we can 
apply the previous lemma with f(x) = xC'(x) and H(z) = e B ( z ' w \ where w is considered as a 
parameter. Then we have 

d B '(x.w)\ ,. [x"]C"(a:,w) ^ _ t [rr"]C fc (x) fc _ ir fe 

ox 1 n^oo [a; n Ja;6'(a;) n-»oo ^-^ [x"JG'(x) f— ' 

and the result follows. 

For the second assertion let us note that B*(x, 1) = B'(x). If we recall that pC'(p) — t and 
tB"{t) = 1, then 

p(l) = pe B ' (T) B"( T ) = P C'{p)t- 1 = 1. 
In order to get an aysmptotic expansion for e^. we have to compute p(u>) explicitly: 

. . t(2D(t) - t)(2D'(t) - l)w 2 ( t(2D{t)~t)w 2 
p{w) — p - — — — — exp tw ' 



2(1 - (2D(t) - t)w) 2 1 V 2(1 - (2D(t) - t)w) , 

This is a function that is admissible in the sense of Hayman [5T| . Hence, it follows that 



d 



k 



where is given by the equation rkp' (rk) /p(rk) = k and b(w) = w 2 p"(w)/p(w) + wp 1 (w) / p{w) — 
(wp 1 (w) / p(w)) 2 . A standard calculation gives the asymptotic expansion for the coefficients dk- □ 

With the help of the explicit expression for p{w) we obtain the values for small k shown in 
Table [U 



4 Series-parallel graphs 

In this section C(x) and B(x) now denote the GFs of connected and 2-connected, respectively, 
series-parallel graphs. First we recall the necessary results from [7] . The radius of convergence of 
B(x) is R sa 0.1280038. The radius of convergence of C(x) is, as for outerplanar graphs, p = ip(r), 
where ij)(u) = ue~ B W, and r is the unique positive root of ip'{u) = 0. Again we have that "0 is 
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the functional inverse of xC'(x), so that r = pC'(p), and r satisfies tB"(t) = 0. The approximate 
values arc r w 0.1279695 and p w 0.1102133. 

In order to study 2-connected series-parallel graphs, we need to consider series-parallel net- 
works, as in [7J. We recall that a network is a graph with two distinguished vertices, called 
poles, such that the graph obtained by adding an edge between the two poles is 2-connected. 
Let D(x, y, w) be the exponential GF of series-parallel networks, where x, y, w mark, respectively, 
vertices, edges, and the degree of the first pole. Define S(x, y, w) analogously for series networks. 
Then we have 

D(x,y,w) = (1 + yw)e s{x > y ' w) - 1 

S(x,y,w) = (D(x,y,w) - S(x,y,w))xD(x,y, 1), 

The first equation reflects the fact that a network is a parallel composition of series networks, and 
the second one the fact that a series network is obtained by connecting a non-series network with 
an arbitrary network (see [3U] for details); the factor D(x, y, 1) appears because we only keep track 
of the degree of the first pole. 



Remark. For the results of the present section, we do not need to take into account the number 
of edges and we could set y = 1 everywhere. However, in the case of planar graphs we do need the 
GF according to all three variables and it is convenient to present already here the full development. 
In the proof of the main result of this section, Theorem 14. b\ we just set y = 1. 

Set E{x, y) = D(x, y, 1), the GF for series-parallel networks without marking the degree of the 
root, which satisfies (see [7J) the equation 

( 1 + E(x,y) \ = xE(x,yf 

V 1 + 2/ J l + xE(x,y)- W 
From the previous equations it follows that 

/ 1 + D(x,y,w)\ xE(x,y)D(x,y,w) 

lo s — n = — rz — wi — \ — • ( 10 ) 

\ 1 + yw J l + xE(x,y) 

Let now B'(x,y) be the GF for 2-connected series-parallel graphs, where the root bears no 
label and has degree k, and where y marks edges. Then we have the following relation. 

Lemma 4.1. 

dB'{x,y,w) s(x,y,w) 

w = xywe y ; . 

ow 

Proof. We have wdB'(x,y,w)/dw = Y2k>i kB*(x,y)w k . The last summation enumerates rooted 
2-connected graphs with a distinguished edge incident to the root, and of these there as many as 
networks containing the edge between the poles (this corresponds to the term e s (x,y,™)y T ne degree 
of the root in a 2-connected graph corresponds to the degree of the first pole in the corresponding 
network, hence the equation follows. □ 

From the previous equation it follows that 

».*->-.»/«*"'** <»> 

Our next task is to get rid of the integral and to express B* in terms of D. Recall that E(x, y) = 
D(x,y,l). 

Lemma 4.2. The generating function of rooted 2-connected series-parallel graphs is equal to 

u . ( x ( ' n , n xE(x,y) ( D(x,y,w) 
B (x, y,w)=x \D(x, y, w) - Y+^Ejxy) ^ ' Vl W ' I 2 
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Proof. We use the techniques developed in [THJ \7\ in order to integrate (fTTj) in closed-form. 



du; = / = y log(l + yw) + / dw 

J l + yw J 1 + yw 

Now we integrate by parts and 

/ ^— dw = y- 1 log(l + yw)L> - / y -1 log(l + yw)—dw. 

J l + yw J dw 

For the last integral we change variables t — D(x,y,w) and use the fact that log(l + yw) = 
log(l + t) - xEt/{l + xE). We obtain 

r qd r D xE f D 

Jlo g (l + y W )—d W = J o logCl + t)*-^-^^ tdt. 

Now everything can be integrated in closed form and, after a simple manipulation, we obtain the 
result as claimed. □ 

In order to prove the main results in this section we need the singular expansions of D(x,y) 
and B(x,y) , for a fixed value of y, near the dominant singularity R(y). 

Lemma 4.3. For \w\ < 1 and for fixed y (sufficiently close to 1) the dominant singularity of the 
functions E(x,y), D(x,y,w), and B* (x,y,w) (considered as functions in x) is given by x — R(y), 
where R(y) is an analytic function in y with R = i?(l) ~ 0.1280038. Furthermore, we have the 
following local expansion: 

E(x, y) = E (y) + E 1 (y )X + E 2 (y)X 2 + ■ ■ • , 
D(x,y,w) = D {y,w) + D 1 (y,w)X + D 2 {y,w)X 2 + ■ ■ ■ , 
B'(x,y,w) = B {y,w) + B 1 (y,w)X + B 2 {y,w)X 2 + ■ ■ ■ , 



where X = yfl — x/R(y). 

The functions R(y), Ej(y), Dj(y, w), and Bj(y, w) are analytic in y resp. in w and satisfy the 
relations 

E (yf = ( log l + im_ Eo{y) 



E (y)-1 V. i + R(v) 

y/1 - 1/E (y) - 1 



Riy)= E 0(y) 
Ei(y) = - 



2R(y)E (y) 2 (l + R(y)E (y)) 2 



(2R(y)E Q (y) + R(y) 2 E (y)^ + 2i?(y)(l + R(y)E (y)) 



Do(y, w) = (1 + yw) exp , D {y, w 



Di(y,w) 
B (y,w) 
Bi (y,w) 



l + R(y)E (y) 
D (y,w)E 1 (y)R(y)(D (y,w) + 1) 
(R(y)E (y)D (y,w)-l)(l + R(y)E Q (y)y 
R(y)D (y,w)(R(y)E (y)D (y,w) -2) 
2(1 + R(y)E (y)) 

2(1 + R(y)E (y)) 2 ' 
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Proof. Since E(x,y) satisfies Equation @ it follows that the dominant singularity of E(x,y) is 
of square-root type and there is an expansion of the form E(x,y) = Eo(y) + E\{y)X + 0(X 2 ), 
with X = yl — x/R(y), and where R(y) and Ej(y) are analytic in y; compare with [TJ 111] . 
Furthermore, if wc set 

^^2 



$0, V, z) = (1 + y) cxp — ) -z-1 

\ 1 + XZ J 

then and Eo(y) satisfy the two equations 

$(R(y),y,E (y))=0 and $ z (R(y), y, E (y)) 
and Ei(y) is then given by 

f2R(y)<S> x (R(y),y,E (y))^ 1/2 



Ei(y) = - 



Next observe that for \w\ < 1 the radius of convergence of the function x i— > D(x,y,w) is 
surely > |i?(y)|. However, D(x,y,w) satisfies Equation (| 10|) . which implies that the dominant 
singularity of E{x,y) carries over to that of D(x,y,w). Thus, the mapping x i— > D(x,y,w) has 
dominant singularity R(y) and it also follows that D(x, y, w) has a singular expansion of the form 
D(x,y,w) = Do(y,w) + Di(y,w)X + 0(X 2 ). Hence, by Lemma [4.11 we also get an expansion 
for B'(x,y,w) of that form. Finally the relations for Dq : D\ and Bq,Bi follow by comparing 
coefficients in the corresponding expansions. □ 

Theorem 4.4. Let dk be the limit probability that a vertex of a two-connected series-parallel graph 
has degree k. then 

i \ a k Bx(l,w) 
v(w) — > di.w = -. 

k>l v ' 

Obviously, p(l) = 1, so that the dk are indeed a probability distribution. We have asymptoti- 
cally, as k — > oo, 

dk-c-k-^W 
where c ~ 3.7340799 is a computable constant and 

g= ((1 + 1/^(1)^0(1))) e -V(i+*(i)Mi)) _ „ 0.7620402. 
Proof. First observe that 

[x n ]B'(x,l,w) 
p(w) = hm - — ^ ; f. 

Fy ' n->oc [x n ]B'( X ,l,l) 

However, from the local expansion of B'(x, 1, w) that is given in Lemma |4~31 fand by the fact that 
B'(x,l,w) can be analytically continued to a A-region; see Section^ it follows that 

[x-]B'{x, l,w) = -^S^ln-^R(1)- (l + O (£)) • 

Hence, p(w) = B 1 (l,w)/B 1 {l,l). 

Next observe that Lemma 14.31 provides Bi(l,w) only for \w\ < 1. However, it is easy to 
continue B\{\,w) analytically to a larger region and it is also possible to determine the dominant 
singularity of B\{1, w), from which we deduce an asymptotic relation for the coefficients of p(w) = 
B 1 (l,w)/B 1 (l,l). 

For this purpose first observe from Lemma 14.31 that D^{y,w) satisfies a functional equation 
which provides an analytic continuation of the mapping w t— » Dv(y,w) to a region including 
the unit disc. In addition, it follows that there exists a dominant singularity wo(y) and a local 
expansion of the form 

D (y,w) =D 00 {y) + D m {y)W + D a2 {y)W 2 + ■ ■ ■ , 
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9(jf,w,z) = (l + yw)exp[ R( - y } E ^ y ] z \ 



where W = y/T- w/wo(y). Furthermore, if we set 

l + R(y)E (y) 
then wo(y) and Doo(y) satisfy the equations 

V(y,w (y),Doo(y))=0 and ^ z (y,w {y),D 00 {y)) = 0. 

Hence 

and " o(y) = K 1 + 55W) exp (" i + R(y)Mv)) ~ ? 



Finally, with the help of Lemma [4.31 it also follows that this local representation of Do(y,w) 
provides similar local representations for D\, Bo, and B\\ 

D 1 {y,w) = + D 10 (y) + D n (y)W + ■ ■ ■ , 

B (y,w) = B 00 (y) + B 02 (y)W 2 + B 03 (y)W 3 + ■ ■ ■ , 
Bi(y, w) = B w (y) + B n (y)W + B 12 (y)W 2 + • ■ ■ , 



where W = yl — w/wo{y) is as above. Hence, all functions of interest Do, D\, Bo, B\ can be 
analytically continued to a A-rcgion, and the asymptotic relation for dk follows immediately. 
Since i«o(l) is the dominant singularity, we have q = l/wo(l). □ 

The next theorem provides the degree distribution in series-parallel graphs. This result has 
been obtained independently in [6], and again our respective results agree. 

Theorem 4.5. Let dk be the limit probability that a vertex of a connected series-parallel graph 
has degree k. then 

I \ T fc ® B'(x,l,w) I 

p{w) = 2_^d k w =p-—e ( ' > \ x =pc'(p) , 
fc>i 

where B* is given by Lemma \4.2\ and Equations UU\) and 0). 

Moreover p(l) = 1, so that the dk are indeed a probability distribution. We have asymptotically, 
as k — > oo ; 

d k ~c-k- 3 / 2 q k , 
where c ks 3.5952391 is a computable constant and 

q=({l + 1/(tE(t, 1))) e -VrS(r,i) _ A' 1 „ o J50416L 

Proof. The proof of the first statement is exactly the same as for Theorem 13.31 Again, we know 
that pC'(p) — t f» 0.127 is larger than the radius of convergence p w 0.110 of C(x), so that 
Lemma |3~T1 applies. The proof that p(l) = 1 is also the same. 

Recall that r < R(l). Hence the dominant singularity x = R(l) of the mapping x i— ► B'(x, 1, w) 
will have no influence to the analysis of p(w). Nevertheless, since 

e B'(x,l,w) = c B'(x,i,w) dB'(x,l,w) 

dx dx 

we have to get some information on D(x, 1, w) and its derivative dD(x, 1, w)/dx with x = r. 

Let us start with the analysis of the mapping w > D{t,\,w). Since D{x,y,w) satisfies 
Equation (fTT)|) it follows that D(t,1, w) satisfies 



D(t, 1, w) = (1 + w) exp 



f tE{t,1)D{t,1,w) 
\ 1 + tE(t, 1) 
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Hence there exists a dominant singularity Wi and a singular expansion of the form 

D(t, 1, to) = D + DiW + D 2 W 2 + ■ ■ ■ , 



where W = \J\ — w/ w\ . Furthermore, if we set 

tE(t, 1)z 



^(w, z) = (1 + w) exp 



1 + tE{t,1] 



then wi and Dq satisfy the equations 

S(wi,£) o ) = and E z (w 1 ,D )=0. 

Consequently, 

Dq = and w\ = I 1 H - r- I exp 



rE(r, 1) V tE(t, I) J V 1 + tE{t, 1) 

Next, by taking derivatives with respect to x in (|10p we obtain the relation 

dD(x,l,w) _ (1 + D(x, 1, w))D(x, 1, w)(E(x, 1) +n£ x (a;,l) 
3x ~~ (x.EO, l)£>(x, 1, tu) - 1)(1 + xE(x, 1)) ' 

Thus if we set x = r and insert the singular representation of D (r, 1 , w) , it follows that dD (^> w ) \ x=r 
has a corresponding singular representation too. By Lemma 14.21 we get the same property for 
dB ^dx ^ \x=t and finally for 

P ■ ^ e B '<-*>W | X=T - Co + C X W + C 2 W^ 2 + • • • . 

This implies the asymptotic relation for dk with q = 1 /wi . □ 

In this case, we obtain an expression for p(w) in terms of the functions E(x, 1) and D(x, 1, w) 
and their derivatives. The derivatives can be computed using Equations ^ and (jlOp as in the 
previous proof. Expanding p(w) in powers of w we obtain the approximate values for small k 
shown in Table [T] 

Remark. We have d\ = p. Also, there is an easy relation between d\ and d 2 , namely 

d 2 = di(2n), 

where nn is asymptotically the expected number of edges in series-parallel graphs. This is shown 
in [551 Thm. 4.10] for planar graphs, phrased in terms of the average degree; the only property 
required is that subdividing an edge preserves planarity, which is also true in the case of series- 
parallel graphs (but not for outcrplanar graphs). The value of k « 1.61673 was determined in [7] 
and one can check that the relation holds. 



5 Quadrangulations and 3-connected planar graphs 

From now on and for the rest of the paper, all generating functions are associated to planar graphs. 
The goal of this section is to find the generating function of 3-connected planar graphs according 
to the degree of the root. This is an essential ingredient in the next section. 

First we work out the problem for simple quadrangulations, which are in bijection with 3- 
connected maps. In order to do that we must revisit the classical work of Brown and Tutte [5] on 
2-connected (non-separable) maps. Finally, using the fact that a 3-connected planar graph has a 
unique embedding in the sphere, we finish the job. 
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5.1 Simple quadrangulations 

A rooted quadrangulation is a planar map where every face is a quadrangle, and with a distinguished 
directed edge of the external face, which is called the root edge of the quadrangulation. The root 
vertex of the quadrangulation is the tail of the root edge. A diagonal is an internal path of 
length 2 joining two opposite vertices of the external face. A quadrangulation is simple if it has 
no diagonal, every cycle of length 4 other than the external one defines a face, and it is not the 
trivial map reduced to a single quadrangle. In Section 5 of [26| it is shown how to count simple 
quadrangulations. Here we extend this result to count them also according to the degree of the 
root vertex. 

A quadrangulation is bipartite and connected, so if we fix the colour of the root vertex there is 
a unique way of 2-colouring the vertices. We call the two colours black and white, and we assume 
that the root is black. Diagonals arc called black or white according to the colour of the external 
vertices they join. 

Let F(x,y,w) be the GF of rooted quadrangulations, where the variables x, y and w mark, 
respectively, the number of black vertices minus one, the number of white vertices minus one, and 
the degree of the root vertex minus one. Generating functions for maps are always ordinary, since 
maps are unlabelled objects. 

The generating functions Fn, Fb and Fw are associated, respectively, to quadrangulations 
with no diagonal, to those with at least one black diagonal (at the root vertex), and to those with 
at least one white diagonal (not at the root vertex). By planarity only one of the two kinds of 
diagonals can appear in a quadrangulation; it follows that 

F(x,y,w) = F N (x,y,w) + F B (x,y,w) + F w (x,y,w). 

A quadrangulation with a diagonal can be decomposed into two quadrangulations, by considering 
the maps to the left and to the right of this diagonal. This gives rise to the equations 

F (x , y , iv) 

F B (x,y,w) = {F N {x,y,w) + F w (x,y,w)) 



F w (x,y,w) = (F N (x,y,w) + F B (x,y,w)) 



x 

F(x,y,l) 



II 



In the second case, only one of the two quadrangulations contribute to the degree of the root 
vertex; this is the reason why the term F(x,y, 1) appears. The x and the y in the denominators 
appear because the three vertices of the diagonal are common to the two quadrangulations. Since 
we are considering vertices minus one, we only need to correct the colour that appears twice at 
the diagonal. Incidentally, no term w appears in the equations for the same reason. 

Let us write F = F(x,y,w) and F{\) = F(x,y,l). From the previous equations we deduce 
that 

F 

x 

F(l) 



F = Fn + Fb+F w = (Fn + F b ) - 
F = Fn + Fb + Fw = (Fn + Fw) yL 
so that 

F + F N = (F N + F B ) + (F N + F w ) = F 



II 



1 + F/x l + F(l)/y 
and finally 

F N =F ( l —— + 1 —--t- - 1 I • (12) 

\l + F/x l + F(l)/y J V ' 

Now we proceed to count simple quadrangulations. We use the following combinatorial de- 
composition of quadrangulations with no diagonals in terms of simple quadrangulations: all quad- 
rangulations with no diagonals, with the only exception of the trivial one, can be decomposed 
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uniquely into a simple quadrangulation q and as many quadrangulations as internal faces q has 
(replace every internal face of q by its corresponding quadrangulation). 
Let 

Q(x, y,w) = ( H.j,kX l y :i w k 

be the GF of simple quadrangulations, where x, y and w have the same meaning as for F. We 
notice that this GF is called Q* N in [26] . We translate the combinatorial decomposition of simple 
quadrangulations into generating functions as follows. 

k / F(l)^ i+j ~ 1 ~ k 
xy/ V xy 



/ F\ ( F(l) 
F N (x, y, w) - xyw = q_i^,kX l y 3 ( — ) ( 

= qi <j> k 



xy {F(1)Y fF(l)\ j ( F 



m v y J v * j \m 
*y o^£W i n i (13) 



where we are using the fact that a quadrangulation counted by qi.j : k has i+j + 2 vertices, — 1 
internal faces, and k of them are incident to the root vertex. 

At this point we change variables as X = F(l)/y, Y = F(l)/x and W = F/F(l). Then 
Equations (11211 and (|T3j can be rewritten as 



V Q(X,Y,W) = F N -xyw = F( - \ / + — 1 -l)-xyw, 



F(l)" v ' " \l + F/x • l + F{l)/y 

Q(X,Y,W) = XYW^ r ^ W¥ + T ^-l^-F(l)w. (14) 

The last equation would be an explicit expression of Q in terms of X, Y, W if it were not for the 
term F(l)w = F(x, y, l)w. In [55] it is shown that 



RS 



F (l) = -i , 5 I 5^3 . ( 15 ) 



where R = R(X, Y) and S(X, Y) are algebraic functions defined by 

R = X(S+1) 2 , S = Y(R+1) 2 . (16) 

Hence it remains only to obtain an expression for w = w(X,Y,W) in order to obtain an explicit 
expression for Q. This is done in the next subsection. 



5.2 Rooted non-separable planar maps 

In [9] the authors studied the generating function h(x,y,w) of rooted non-separable planar maps 
where x, y and w count, respectively the number of vertices minus one, the number of faces minus 
one, and the valency (number of edges) of the external face. We notice that the variable z is used 
instead of w in [9]. There is a bijection between rooted quadrangulations and non-separable rooted 
planar maps: black and white vertices in the quadrangulation correspond, respectively, to faces 
and vertices of the map; quadrangles become edges; and the root vertex becomes the external face, 
and its degree becomes its valency. As a consequence h(y,x,w) = wF(x,y,w), where the extra 
factor w appears because in F we are counting the degree of the root vertex minus one. It follows 
that Equation (3.9) from [3] becomes 

(1 - io)(l - yw)wF = -w 2 F 2 + (-xw + wF(l))wF + xw 2 (x(l - w) + F(l)). 



14 



By dividing both sides by F(l) 2 and rewriting in terms of X = F(l)/y, Y = F(i)/x and W 
F/F(l), we obtain 

, / 1 w\ oof l\o w 2 ( 1 

(1 - w) \ — — - — )wW= -~w 2 W 2 





w 2 ( 


(1-7) 









F{1) X, 

(1 - w)(X - wF(l))wW _ w 2 (-XYW 2 + XYW - XW + 1 - w + X 
XF(1) ~ XY ' 

Y(l - w)(X - wF(l))W = wF{l){-XYW 2 + XYW - XW +1-W + X). (17) 

Observe that this is a quadratic equation in w. Solving for w in (|17p and using (|15p and (|16p we 
get (the plus sign is because T* has positive coefficients in coming Theorem 15 - 1[) 



Wl (R, S,W) + (R~W + 1) y/w 2 (R, S, W) 



2(5 +l) 2 (SW + R 2 + 2R+1 
where Wi(R, S, W) and W2(R, S, W) are polynomials given by 



iui = - RSW 2 + W{1 + 4S + 3RS 2 + 5S 2 + R 2 + 2R + 2S 3 + 3R 2 S + 7RS) (19) 

+ {R+l) 2 (R + 2S+l + S 2 ), 
w 2 =R 2 S 2 W 2 - 2WRS(2R 2 S + 6RS + 2S 3 + 3RS 2 + 5S 2 + R 2 + 2R + AS + 1) (20) 

+ {R+1) 2 (R + 2S+1 + S 2 ) 2 . 

The reason we choose to write w as a function of (R, S, W) instead of (X, Y, W) will become clear 
later on. 

Thus, together with Equations (| 14[) and (|15[) . we have finally obtained an explicit expression 
for the generating function Q(X,Y,W) of simple quadr angulations in terms of W and algebraic 
functions R(X,Y) and S(X,Y). 



5.3 3-connected planar graphs 

Let T*(x, z, w) be the GF of 3-connected planar graphs, where one edge is taken as the root and 
given a direction, and where x counts vertices, z counts edges, and w counts the degree of the tail 
of the root edge. Now we relate T* to the GF Q(X, Y, W) of simple quadrangulations. 

By the bijection between simple quadrangulations and 3-connected planar maps, and using 
Euler's relation, the GF xwQ{xz, z,w) counts rooted 3-connected planar maps, where z marks 
edges (we have added an extra term w to correct the 'minus one' in the definition of Q). 

According to Whitney's theorem 3-connected planar graphs have a unique embedding in the 
sphere. As noticed in [3], the are two ways of rooting an embedding of a directed edge- rooted 
graph in order to get a rooted map, since there are two ways of choosing the root face adjacent to 
the root edge. It follows that 

XVJ 

T'(x,z,w) = —— Q(xz,z,w). (21) 

Theorem 5.1. The generating function of directed edge-rooted 3-connected planar graphs, where 
X,Z,W mark, respectively, vertices, edges, and the degree of the root vertex, is equal to 



^ # x 2 z 2 w 2 I 



(u + l) 2 W\{U, V,w) + (U — W + l)y / W2{u, V, w)^j 



1 + ivz 1 + xz 2w(vw + u 2 + 2u + 1)(1 + u + v) 3 

(22) 

where u and v are algebraic functions defined by 

u = xz(l + v) 2 , v = z{l + u) 2 , (23) 
and W\(u, v, w) and 11)2(11, v, w) are given by jT9j) and \20\) replacing R, S, W by u, v, w, respectively. 
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Proof. Combine Equation lj2"T]). together with Equations flij). P^|) . and ([H]). □ 

When we set w = 1 in Equation (|2"2")l we recover the GF of edge-rooted 3-connected planar 
graphs without taking into account the degree of the root vertex. 

6 Planar graphs 

This section is divided into three parts. First we obtain an explicit expression for B'(x,y,w), 
the generating function of rooted 2-connected planar graphs taking into account the degree of the 
root. Secondly, we compute singular expansions at dominant singularities for several generating 
functions. And finally we obtain the asymptotic degree distribution in random planar graphs. 

6.1 2-connected planar graphs 

Let B*(x,y,w), the generating function of rooted 2-connected planar graphs taking into account 
the degree of the root. As for series series-parallel graphs we have to work with networks. 

Let T'{x, z, w) be the GF for directed edge-rooted 3-connected planar maps as in the previous 
section. As in Section 21 we denote by D(x, y, w) and S(x, y, w), respectively, the GFs of (planar) 
networks and series networks, with the same meaning for the variables x, y and w. 

Lemma 6.1. We have 

D(x,y,w) 



D(x,y,w) = (l + yw)exv[S(x,y,w) + -^ F - -T" I x,E{x,y), 

\ x z D(x,y,w) \ E{x,y) 

S{x,y,w) = xE(x,y)(D(x,y,w) - S(x,y,w)) , 
where E(x,y) = D(x,y 7 1) is the GF for planar networks (without marking the degree of the root). 

Proof. The proof is a variant of the equations developed by Walsh [3D], taking into account the 
degree of the first pole in a network. The main point is the substitution of variables in TV an 
edge is substituted by an ordinary planar network (this accounts for the term E{x,y)), except if 
it is incident with the first pole, in which case it is substituted by a planar network marking the 
degree, hence the term D(x, y, w) (it is divided by E(x, y) in order avoid overcounting of ordinary 
edges). □ 



As in Lemma l4.1| and for the same reason, we have 

k>l 

(at \ i 1 rr. ( jpf n D(x,y,w) 

= xywexp[b(x,y,w) + - T =- rT \x,E{x,y), —- r- 

V x z D(x,y,w) \ E{x,y) 

Lemma 6.2. The generating function of rooted 2-connected planar graphs is egual to 

( xED ( D\\ 1 + D . . 1 f D T m (x,E,t/E) 

where for simplicity we let D — D(x, y, w) and E = E{x, y). 
Proof. We start as in the proof of Lemma 14.21 

/ } 7- — dw = jT 1 log(l + yw) + 2T 1 log(l + yw)D - f y" 1 log(l + yw)^- dw. 
J 1 + yw J ow 

For the last integral we change variables t = D(x, y, w) and use the fact that 
log(l + yw) = log(l + D) - { ^ E) -^T'(x, E, D/E). 



dt, (24) 
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Wc obtain 



, , f° , „ x, xE f D , 1 f D T*(x,E,t/E) 1 

log(l + iH 3^ ^ - y Q Iog(l + 1) dt - ^ tdt + -^J t 



On the right-hand side, all the integrals except the last one are elementary. Now we use 

and after a simple manipulations the result follows. □ 



B'(x,y,w) =xy J y 



In order to get a full expression for B*(x,y,w), it remains to compute the integral in the 
formula of the previous lemma. 

Lemma 6.3. Let T'(x, z, w) be the GF of 3-connected planar graphs as before. Then 

T'(x, z, t) _ x 2 (z 3 xw 2 - 2wz - 2xz 2 w + (2 + 2xz) log(l + wz)) 
t - 4(1 + xz) 

uvx ( w(2u 3 + {6v + 6)u 2 + (6v 2 - vw + lAv + 6)u + Av 3 + 10w 2 + 8v + 2) 



2(1 + u + v) 3 \ Av(v + I) 2 

(1 + u)(l + u + 2v + v 2 )(2u 3 + (4i> + 5)u 2 + (3v 2 + 8v + 4)u + 2v 3 + 5v 2 + Av + 1) 

Auv 2 (v + l) 2 

~Q(2u 3 + (4v + 5)u 2 + (3v 2 - vw + 8v + A)u + 5v 2 + 2v 3 + 4v + 1) 



4uv 2 (v + I) 2 

(l + M) 2 (l + M + i;) 3 log(Qi) 
+ 2v 2 (l + v) 2 

(u 3 + 2u 2 + u - 2v 3 - 4u 2 - 2v)(l +u + v) 3 log(Q 2 ) 
+ 2v 2 (l + v) 2 u 

where the expressions Q, Qi and Q2 are given by 



= u 2 v 2 w 2 - 2uvw{u 2 (2v + 1) + u{3v 2 + 6v + 2) + 2v 3 + 5v 2 + Av + 1) 
+ (l + u) 2 (u + {v + l) 2 ) 2 

1 



(-uvw(u 2 + u(v + 2) + 2v 2 + 3v + 1) 



2(wv + {u + l) 2 ) 2 (v + l)(u 2 + u{v + 2) + (v + l) 2 ) 
+(u + l)(u + v + l)y/Q + (u + l) 2 (2u 2 {v + 1) + u{v 2 + 3v + 2) + v 3 + 3v 2 + 3v + lj 
-wuv + u 2 (2v - 1) + u(3v 2 + 6v + 2) + 2v 3 + 5v 2 + Av + 1 - VQ 



2v(u 2 + u(v + 2) + (v + l) 2 ) 



Proof. We use Equation l]22[) to integrate T'(x,z,w)/w. Notice that neither u nor v have any 
dependence on w. We have used Maple to obtain a primitive, to which we have added the 
appropriate constant c(x, z) to ensure that the resulting expression evaluates to when w = 0. 

A key point in the previous derivation is that, by expressing w(X, Y, W) (sec Equation (|17[) ) in 
terms of R, S instead of X, Y (see Equation (fI5|)). we obtain a quadratic polynomial w 2 (R, S, W) 
in terms of W inside the square root of T'(x, z, w) in Equation (|22p . Otherwise, we would have 
obtained a cubic polynomial inside the square root, and the integration would have been much 
harder. □ 

Combining Lemmas 16.21 and 16.31 we can produce an explicit (although quite long) expression 
for B'{x,y,w) in terms of D(x, y, w), E(x,y), and the algebraic functions u(x,y),v(x,y). This 
is needed in the next section for computing the singular expansion of B'(x,y,w) at its dominant 
singularity. 
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6.2 Singular expansions 

In this section we find singular expansions of T*(x, z, w), D(x, y, w) and B*(x, y,w) at their dom- 
inant singularities. As we show here, these singularities do not depend on w and were found in [3] 
and [19 . But the coefficients of the singular expansions do depend on w, and we need to compute 
them exactly in each case. In the next section we need the singular expansion for B* , but to 
compute it we first need the singular expansions of u, v, T* and D (for u and v, see also [UG2). 

Lemma 6.4. Let u = u(x, z) and v = v(x, z) be the solutions of the system of equations u = 
xz{l + v) 2 and v = z(l + u) 2 . Let r(z) be given explicitly by 

_ Up(z) 

( >~ z(l + z(l + u (z)) 2 ) 2 ' [Zb} 

where 



w °(*) = -3 + V9 + 3? 
Furthermore, let t(x) be the inverse function of r(z) and let uq(x) = u~o(i~(x)) which is also the 
solution of the equation 

(l + u)(3u- l) 3 

x = 7. ■ 

lou 

Then, for x sufficiently close to the positive real axis the function u(x, z) and v(x, z) have a 
dominant singularity at z — t(x) and have local expansions of the form 

u(x, z) = uq(x) + u\{x)Z + u 2 (x)Z 2 + u^(x)Z 3 + 0{Z A ), 
v(x, z) = v (x) + Vl {x)Z + v 2 {x)Z 2 + v 3 (x)Z 3 + 0(Z 4 ), 



where Z = y/l — z/t(x). The functions Uj(x) and Vj(x) are also analytic and can be given 
explicitly in terms of u — Uq(x). In particular we have 



uo(x) = u vo(x) 



1 



3u-l 



, x / — ; r , * 2J2u(u + 1) 

U 1 (x)=-y/2u(u+l) Vl{ X ) = V 1 ' 

Su — 1 

, . (l + u)(7« + l) , , 2m(3 + 5u) 

M x ) = on , o ^ V ^ X > ~ 



2(1 + 3u) K ' (3u- l)(l + 3u)' 

_ (1 + m)(67m 2 + 50m+11)m _ ygu(l + u)(79u 2 + 42u + 7) 

~ 4(l + 3u) 2 V2u 2 + 2u ~ ~ 4(l + 3u) 2 {3u-l)yfuTT+u)' 

Similarly, for z sufficiently close to the real axis the functions u = u(x, z) and v = v(x, z) have 
a dominant singularity x = r{z) and there is also a local expansion of the form 

u(x, z) = u (z) + ui(z)X + u 2 {z)X 2 + 0(X 3 ), 
v(x, z) = v Q (z) + vi{z)X + v 2 {z)X 2 + 0(X 3 ), 



where now X = yl — x/r(z). The functions Uj(z) and Vj(z) are analytic and can be given 
explicitly in terms ofu = uq(z). In particular, we have 



u (z) = u v (z) 



1 



32-1 



„ . , 2ua/1 + u „ , , Auy/T 

u n z ) = 7= v i( z ) 



\/l + 3u (3u- l)Vl + 3u 

„ . , 2(1 + u)u(2u + 1) „ , , 4u(5u 2 + 4u+l) 

M z ) = n : o~N2 w) ~ 



(l + 3u) 2 v ' (3u-l)(l + 3u) 



2 



_ , , 2S(10m 3 + 112 2 + 55+1)\/TTS 4m(2u + 1)(11m 2 + 5u + l)vTT~ 



(l + 3u) 7 /2 ^ > (3u-l)(l + 3u) 7 /2 
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Proof. Since u(x,z) satisfies the functional equation u = xz(l + z(l + u) 2 ) 2 , it follows that for 
any fixed real and positive z the function x i— > u(x, z) has a square-root singularity at r(z) that 
satisfies the equations 

$(r(z),z,u) = and $ u (r(z), z, it) = 0, 

where z, u) = u — xz(l + z(l + u) 2 ) 2 . Now a short calculation gives the explicit formula (|25|) 
for r(z). By continuity we obtain the same kind of representation if z is complex but sufficiently 
close to the positive real axis. 

We proceed in the same way if x is fixed and z is considered as the variable. Then t{x), 
the functional inverse of r(z), is the singularity of the mapping z t— > u(x,z). Furthermore, the 
coefficients Ui(x) etc. can be easily calculated. The derivations for v{x,z) arc completely of the 
same kind. □ 

Lemma 6.5. Suppose that x and w are sufficiently close to the positive real axis and that \w\ < 1. 
Then the dominant singularity z = t(x) ofT'(x,Z,w) does not depend on w. The singular 
expansion at t(x) is 

T'(x,z,w) =T (x,w) + T 2 (x,w)Z 2 +T 3 (x,w)Z 3 + 0{Z 4 ), (26) 

where Z = y/1 — z/t(x), and the expressions for the Ti are given in the appendix. 

Proof. Suppose for a moment that all variables x,z 7 w are non- negative real numbers and let us 
look at the expression ([22]) for TV The algebraic functions u(x,z) and v(x,z) are always non- 
negative and, since the factor w in the denominator cancels with a corresponding factor in the 
numerator, the only possible source of singularities are: a) those coming from u and v, or b) the 
vanishing of W2(u,v,w) inside the square root. 

We can discard source b) as follows. For fixed u, v > 0, let W2{w) = W2(u, v, w). We can check 
that 

102(1) = (1 + 2u + u 2 + 2v + v 2 + uv - uv 2 ) 2 , 

w' 2 (w) = -2uv ((6 - w)uv + 1 + 2u + u 2 + 4v + 5v 2 + 2v 3 + 3uv 2 + 2u 2 v) . 

In particular 102(1) > and w' 2 (w) < for w € [0,1], thus it follows that W2(w) > in 
w G [0, 1]. Hence the singularities come from source a) and do not depend on w. 

Following Lemma IB~4l (see also [4] and [3]), we have that z = t{x) is the radius of convergence 
of u(x, z), as a function of z. Now by using the expansions of u(x, z) and v(x, z) from Lemma l6.4l 
we obtain (l26l) . 

Finally, by continuity all properties arc also valid if x and w are sufficiently close to the real 
axis, thus completing the proof. □ 

Similarly we get an alternate representation expanding in the variable x. 

Lemma 6.6. Suppose that z and w are sufficiently close to the positive real axis and that \w\ < 1. 
Then the dominant singularity x — r(z) ofT*(x,z,w) does not depend on w. The singular 
expansion at r{z) is 

T'{x, z, w) = f (z, w) + T 2 (z, w)X 2 + T 3 (z, w)X 3 + 0(X 4 ), (27) 

where X = y/l — x/r(z). Furthermore we have 

f (z,w) = T (r(z),w), 

T 2 (z, w) = T 2 (r(z), w)H(r(z), z) ~ T , x (r(z), x)r(z), 
T 3 (z,w)=T 3 (r(z),w)H(r(z),z)V 2 , 

where H(x, z) is a non-zero analytic function with Z 2 = H(x, z)X 2 . 
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Proof. We could repeat the proof of Lemma 16.51 However, we present an alternate approach that 
uses the results of Lemma 16.51 and a kind of singularity transfer. 

By applying the Weicrstrass preparation theorem it follows that there is a non-zero analytic 
function with Z 2 = H(x, z)X 2 . Furthermore, by using the representation x — r(z)(l — X 2 ) and 
Taylor expansion we have 

H(x, z) = H(r(z), z) - H x (r(z),z)r(z)X 2 + 0{X 4 ), 
Tj{x, w) = Tj(r(z), w) - T hX {r{z),w)r{z)X 2 + 0(X 4 ). 

Hence, Lemma (|6.5p directly gives proves the result. In fact, H can be computed explicitly and is 
equal to H(x,t(x)) = (1 + 3u)/2u. □ 

The next result is Proposition 6.3 from [12j . and is needed in order to guarantee that singular 
expansions of the desired kind exist. 

Theorem 6.7. Suppose that F(x,y,u) has a local representation of the form 

Fix, y, u) = g(x, y, u) + h(x, y,u)[l r ) (28) 

with functions g(x,y,u), h(x,y,u), r(x,u) that are analytic around (xo,yo,uo) and satisfy 
g y (x ,y ,u Q ) ^ 1, h(x ,y ,u ) ^ 0, r(x ,u )^Q, r x (x , u ) ^ g x (x , yo, u ). 
Furthermore, suppose that y = y(x, u) is a solution of the functional equation 

y = F(x,y,u) 

with y(xQ,uo) = yo- Then y[x,u) has a local representation of the form 

( * \ 3/2 

y{x,u) = gi(x,u) + hi{x,u) I 1 - -^-y I , (29) 

where gi(x, u), h\(x, u) and p(u) are analytic at (xq, uq) and satisfy h\(xo, uo) ^ and p(uq) — xq. 

Lemma 6.8. Suppose that y and w are sufficiently close to the positive real axis and that \w\ < 1. 
Then the dominant singularity x = R(y) of D{x,y,w) does not depend on w. The singular 
expansion at R(y) is 

D(x, y, w) = D (y, w) + D 2 (y, w)X 2 + D 3 (y, w)X 3 + 0(X 4 ), (30) 

where X = \J\ — x/R(y), and the expressions for the Di are given in the appendix. 

Proof. We use the system of equations from Lcmma l6.ll If we set w = 1 and substitute S(x, y, 1) = 
xE(x, y) 2 1 (1 + xE(x, y)) then we get a single equation for E(x, y) = D{x, y, 1): 

E { ,, „) = <! + 9 )cxp (y^fgy + J^jnx, »).!))-!■ (3D 

Now we use the singular expansion from Lemma 16.61 and Theorem 16.71 (which is Proposition 6.3 
of [12]) to conclude that E(x, y) has an expansion of the form 

E(x,y) =E (y) + E 2 (y)X 2 + E 3 (y)X 3 + O(X 4 ). (32) 

Finally we reconsider Lemma 16.11 and after substituting S(x,y,w) = xE(x,y)D(x,y,w)/(l + 
xE{x,y)) wc get a corresponding (single) equation for D(x,y,w): 

rv \ ft , \ f xE(x,y)D(x,y,w) 1 ( D(x,y,w)* 

D(x, y, w) = (l + yw) exp , r + zT x, E(x, y), — - 

V 1 + xE{x, y) x z D{x,y,w) \ E(x,y) 

(33) 

Now a second application of Lemma RDjI and Theorem 16.71 yields the result. Note that the singu- 
larity does not depend on w. □ 
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Lemma 6.9. Suppose that y and w are sufficiently close to the positive real axis and that \w\ < 1. 
Then the dominant singularity x = R(y) of B'(x,y,w) does not depend on w, and is the same as 
for D(x,y,w). The singular expansion at R(y) is 

B'(x, y, w) = B (y, w) + B 2 (y, w)X 2 + B 3 (y, w)X 3 + 0{X A ), (34) 

where X = y/l — x/R(y), and the expressions for the Bi are given in the appendix. 

Proof. We just have to use the representation of B*(x,y,w) that is given in Lemma 16.21 and 
Lemma 16.31 and the singular expansion of D(x,y,w) from Lemma 16.81 □ 



6.3 Degree distribution for planar graphs 

We start with the degree distribution in 3-connected graphs, both for edge-rooted and vertex- 
rooted graphs. 

Theorem 6.10. Let dk be the limit probability that a vertex of a three- connected planar graph has 
degree ft, and let be the limit probability that the tail root vertex of an edge-rooted (where the 
edge is oriented) three- connected planar graph has degree k. Then 



E 



e k w 



k _ T 3 (r,w) 



where T 3 (x,w) is given in the appendix and r = r(l) = (7y/7 — 17)/32 is explicitly given by \2_ 
Obviously the are indeed a probability distribution. We have asymptotically, as k — > oo, 



where c w 0.9313492 is a computable constant and q = l/(uo + 1) = \fl — 2, and where uq = 
u{r) = (V7-1J/3. 
Moreover, we have 

d k =a— ~ca- k~ 1/2 q k , 
k 

where 

(3uo-l)(3uo + l)(uo + l) V7 + 7 
a = = — ^ — 

UQ 2 

is the asymptotic value of the expected average degree in 3-connected planar graphs. 

We remark that the degree distribution in 3-connected planar maps counted according to the 
number of edges was obtained in [2|. The asymptotic estimates have the same shape as our dk, 
but the corresponding value of q is equal to 1/2. 

Proof. The proof uses first the singular expansion (|27p . The representation 

k T 3 (r,w) 

2_£kw = — 

h>\ y ' ' 

follows in completely the same way as the proof of Theorem 14.41 using now Lemma 16. 5[ with the 
difference that now the dominant term is the coefficient of Z 3 . 

In order to characterize the dominant singularity of T 3 (l,w) and to determine the singular 
behaviour we observe that the explicit representation for T 3 contains in the denominator a (dom- 
inating) singular term of the form (— w + u + l) 3 / 2 . Hence, it follows that uq + 1 is the dominant 
singularity and we also get the proposed asymptotic relation for efc. 

For the proof of the second part of the statement, let t Ut k be the number of vertex-rooted graphs 
with n vertices and with degree of the root equal to k, and let s n ^k be the analogous quantity 
for edge-rooted graphs. Let also t n = J2k^ n , k an< ^ s ™ = 12k Sn - k - Since a vertex-rooted graph 
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with a root of degree k is counted k times as an orientcd-edgc-rootcd graph, we have s n< k = kt n ^ 
(a similar argument is used in [22j ) . Notice that ek = lbcas nt k/s n an( l dk = lim 

Using the quasi-powers theorems as in [19) . one shows that the expected number of edges /i n 
in 3-connected planar graphs is asymptotically [i n ~ nn, where k = — t'(1)/t(1), and t(x) is as in 
Lemma 16.41 Clearly s n = 2/i n t n /n. Finally, 2fi n /n is asymptotic to the expected average degree 
a = 2k. Summing up, we obtain 

kd k = ae k - 

A simple calculation gives the value of a as claimed. □ 

Theorem 6.11. Let dk be the limit probability that a vertex of a two-connected planar graph has 
degree k. Then 

B 3 (l,w) 
P ^ = IMM? 

where B 3 (y,w) is given in the appendix. 

Obviously, p(l) = 1, so that the dk are indeed a probability distribution and we have asymptot- 
ically, as k — > oo, 

d k ~ cfc-^V, 
where c ~ 3.0826285 is a computable constant and 

1 _f(to-l){to + 6) 



exp ^ — - - 1 w 0.6734506, 

I- to \ 6*2 + 20*0 ' 



and *o = t(X) ~ 0.6263717 is a computable constant given in the appendix. 



Proof. The representation of p(w) follows from (f34|) . 

Now, in order to characterize the dominant singularity of B 3 (l, w) and to determine the singular 
behaviour we first observe that the right hand side of the equation for Do contains a singular term 
of the form (Do(t— 1) + t) 3 ^ 2 that dominates the right hand side. Hence, by applying Theorem 16. 71 
it follows that Do(l,w) has dominant singularity IU3, where Do (1,^3) = t/(l — t) and we have a 
local singular representation of the form 



D (l, w) = Doo + D 02 W 2 + D 03 W 3 



where W = y/T— w/w 3 and D o = */(l — *). The fact that the coefficient of W vanishes is due 
to the shape of the equation satisfied by D (l,w) (see the appendix). 

We now insert this expansion into the representation for D 2 . Observe that 5 has an expansion 
of the form 

S = S 2 W 2 + S3W 3 +■■■ 

with S2 7^ 0. Thus we have y/S = y/S^W + 0(W 2 ) . Furthermore, we get expansions for S^i, 5*2,2, 
52,3, and 52,4- However, we observe that 52,3(1, w 3 ) = whereas S2,]. (1,103) 7^ 0, 52,2(1} 103) 7^ 0, 
and 52,4(1, w 3 ) 7^ 0. Consequently we can represent 1)2(1, w ) as 

D 2 (l,w) = D 2 ^ X L + D 2 .o + D 2 .iW + +D 22 W 2 + ■ ■ ■ , 
W 

where D\ —i 7^ 0. 

In completely the same way it follows that -©3(1, w) has a local expansion of the form 

d 3 (i, w) = 5 3i _ 3 ^- + +5 3) _ 1 ^ + 5 3j0 + 53,1^7 + ••■ , 
w 3 w 

where 1)3,-3 7^ 0, and the coefficient of W~ 2 vanishes identically. 

These type of singular expansions carry over to B 3 (l,w), and we get 

B 3 (l,w) = B^JL + B 3 , + B 3il W + ■ ■ ■ , (35) 
W 
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where B3.-1 7^ 0. We stress the fact that the coefficients of W 3 and W 2 vanish as a consequence 
of non trivial cancellations. Hence, we obtain the proposed asymptotic relation for the dk- □ 

The following is the analogous of Lemma 13.11 The difference now is that we are composing 
two singular expansions and, moreover, they are of type 3/2. 

Lemma 6.12. Let f(x) — Y^ n >o a nX n /n\ denote the exponential generating function of a sequence 
a n of non-negative real numbers and suppose that f{x) has exactly one dominating singularity at 
x = p of the form 

f(x) = f + hX 2 + f 3 X 3 + G(X 4 ), 



where X = yl — xj p, and has an analytic continuation to the region {i£C: |ar| < p + e}\ {x € 
M. : x > p} for some s > 0. Further, let H(x,z,w) denote a function that has a dominant 
singularity at z = f(p) > of the form 

H(x, z, w) = h (x, w) + h 2 (x, w)Z 2 + h 3 (x, w)Z 3 + 0{Z A ), 



where w is considered as a parameter, Z = \J\ — z/f(p), the functions hj(x,w) are analytic in x, 
and H(x, z, w) has an analytic continuation in a suitable region. 
Then the function 

f H (x)=H(x,f(x),w) 
has a power series expansion fnix) = X)n>o b n x n jn\ and the coefficients b n satisfy 

n ^°° a n fa f 3 V Jo J 

Proof. The proof is similar to that of Lemma [8] and is based on composing the singular expansion 
of H(x, z, w) with that of f{x). Indeed, near x = p, and taking into account that f(p) = fo, we 
have 

f 2 X 2 + f 3 X 3 \ . / f 2 X 2 + f 3 X^ 3/2 



fn(x) = h (x,w) + h 2 (x,w) + h 3 (x,w, . 

\ To J V ./( 

Now note that x = p — pX 2 . Thus, if we expand and extract the coefficient of X 3 and applying 
transfer theorems, we have 

— ~ ■fc n - 5/2 p- n 



n\ T(-3/2) 



and 



h 2 (p, w)f 3 ( f 2 \ 3/2 \ _ 5/2 _„ 



n\ r(-3/2) 

so that the result follows. □ 

Theorem 6.13. Let dk be the probability that a vertex of a connected planar graph has degree k. 
Then 



k>l 

+ eB .( 1 , M )-^( 1 ,l) l + yi,l) Ba(1>||>)| 
-0311, 1) 

where Bj(y,w), j = 0,2,3, are given in the appendix. 

Moreover, p(l) = 1, so that the dk are indeed a probability distribution and we have asymptot- 
ically, as k — > 00, 

d k ~ ck~ 1 ' 2 q k , 

where c f=a 3.0175067 is a computable constant and q is as in Theorem lO. 1 1{ 
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Proof. The degree distribution is encoded in the function 



w k = e B'(xC (x),l,w) 



k>l 

where the generating function xC'(x) of connected rooted planar graphs satisfies the equation 

xC / (x)=xe B '^ c '^ 1 l 
From Lemma 16.91 we get the local expansions 

e B'(x,i,w) = e B (i. w ) ^ + b 2 (i jW }x 2 + B 3 (l lW )X 3 + 0{X 4 )) , 

where X = yl — x/R. Thus, we first get an expansion for xC'(x) 

xC'(x) =R - t -X 2 + - RWA \^ l0 X* + 0(X 4 ), 

V ' 1 + B 2 (1,1) (1 + 5 2 (1,1)) 5 / 2 V ; ' 



where X = — xj p and p is the radius of convergence of C'{x) (compare with [19]). Note 
also that R = pe 3 "^-^. Thus, we can apply Lemma IFTT21 with H(x,z,w) = xe 3 '^' 1 ^ and 
f(x) = xC'{x). We have 

Jo — it, h — - , „ , J3 



1 + 52(1,1)' JA (1 + B 2 (1,1)) 5 / 2 
and 

h (p,w) =pe B ° ( - 1 > w \ h 2 (p,w)=pe B ^B 2 (l,w), h 3 (p, w) = pe B ^B 3 (l, w). 
We can express the probability generating function p(w) as 

\x n ]xC*(x,w) 

hm — : — — . 

n-*oc [X n \xC (X) 

Consequently, we have 

h 2 (p,w) h 3 (p,w) ( f 2 \ 3/2 

p{w) = — 



fo h V fo 

= -e B °^- B °^B 2 (l,w) 

+ B,(l,«,)-B„(l,l) 1 + B 2 (l, l) fl (1 , 
B 3 (l,l) 

The singular expansion for £> 2 (1, w) turns out to be of the form 

B 2 (l,w) = B 2 . Q + B 2A W + ■■■ . 

Hence the expansion ([55)1 for B 3 (l, w) gives the leading part in the asymptotic expansion for p(w). 
It follows that p(w) has the same dominant singularity as for 2-connected graphs and we obtain the 
asymptotic estimate for the dk as claimed, with a different multiplicative constant. This concludes 
the proof of the main result. □ 



7 Degree distribution according to edge density 

In this section we show that there exists a computable degree distribution for planar graphs with 
a given edge density or, equivalently, given average degree. 
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In [19] we showed that, for p in the open interval (1,3), the number of planar graphs with n 
vertices and \_pn\ ls asymptotically 

9,1,1^] ~ c(fi)n- i ^(n) n n\, (37) 

where c(p) and 7(//) are computable analytic functions of p. The proof was based on the fact for 
each p G (1,3), there exists a value y = y(fi) such that the generating function G(x,y) captures 
the asymptotic behaviour of g n ,l^n\ ■ The exact equation connecting p and y is 

P'(y) f oo\ 

~ V^-j = Mi (38) 

where p(y) is the radius of convergence G(x, y) as a function of x. More precisely, the idea behind 
the proof is to weight a planar graph with m edges by y m . If g nm denotes the number of planar 
graphs with n vertices and m edges, then the bivariate generating function 

™n ™n 



G{x, y)=Yl 9nm — y m = 9n(y) 



m,n n>l 



can be considered as the generating function of the weighted numbers g n (y) = ^2 m gnm.y m of 
planar graphs of size n. In addition this weighted model induces a modified probability model. 
Instead of the uniform distribution the probability of a planar graph of size n is now given by 
y m /g n {y), where m denotes the number of edges. It follows from the singular expansion [19] 

G(x, y) = G (y) + G 2 (y)Y 2 + G A (y)Y 4 + G 5 (y)Y 5 + Q(Y 6 ), 



where Y = yl — x/p(y), that g n {y) is asymptotically given by 

g n (y)~g(y)n- 7 / 2 p(y)- n n\ 

with g(y) = G${y) /Y{— 5/2). If y = 1 we recover the asymptotic formula for the number of planar 
graphs of size n. However, since g n {y) is a power series in x with coefficients g nm it is possible to 
compute these numbers by a Cauchy integral: 

1 f , x dy 

g nm = 7T- / 9n{y) 



2iri J\y\= r V m+1 ' 

Suppose that m = fm. Then due to the asymptotic structure of g n (y) the essential part of the 
integrand behaves like a power: 

g n {y)y~ rn ~ g{y)n~ 7 l 2 n\ {y»p(y))- n = g{y)n 7 / 2 n\ cxp (-n log{y^p(y))) . 

Hence, the integral can be approximated with the help of a saddle point method, where the saddle 
point equation ^ \og(y' 1 p(y)) = is precisely (|3"8j) . This leads directly to ([3T|) : see [12] for details. 
Informally, one has to append just a saddle point integral at the very end of the calculations. 

We can use exactly the same approach for the degree distribution. We consider the weighted 
number of rooted planar graphs of size n, where the root has degree k. For fixed y, the correspond- 
ing generating function G'{x, y, w) has the property that the radius of convergence of G*(x, y, w) 
does not depend on w and is, thus, given by p(y). The asymptotic equivalent for <?*(?/, w) contains 
the factor p(y) n and consequently it is possible to read off the coefficient of y m with the help of 
a saddle point method as above. Note that we have computed in the appendix the coefficients of 
the singular expansion of B(x, y, w) as a function not only of w but also of y. Hence, Theorem 
16.131 extends directly to the following. 

Theorem 7.1. Let p G (1,3) and let d^^k be the probability that a vertex of a connected planar 
graph with edge density p has degree k. Let y be the unique positive solution of t38]). Then 

Y, d ^w k - -e B °^- B °^B 2 (y,w) 

k>l 

B 3 (y,l) 3ly ' h 

where Bj(y,w), j — 0,2,3 is given in the appendix. 
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The probabilities d^.k can be computed explicitly using the expressions for the B(y, w) given 
in the appendix. As an illustration in Figure [T] we present a cumulative plot for k = 1, . . . , 10. 
Each curve gives the probability that a vertex has degree at most k. The abscissa is the value of 
/i £ ( 1, 3) , while the ordinate gives the probability. The bottom curve corresponds to k = 1 and 
the top curve to k < 10. The vertical line is the value k « 2.21 such that nn is the asymptotic 
expected number of edges in planar graphs [19j ; since the number of edges is strongly concentrated 
around nn, the probabilities in these abscissa correspond to the cumulating values in the third 
line of Table [U namely 0.0367284,0.1993078,0.4347438,0.6215175,0.7510198,0.8372003,... 




Figure 1: Cumulative degree distribution for connected planar graphs with \in edges, /i £ (1,3) 
and k = 1, . . . , 10. 

Figure [2] shows the data for 2-connected and 3-connected planar graphs. Notice that in a 

2- connected graph the degrees are at least 2, in a 3-connected graph they are at least 3, and a 

3- connected graph has at least 3n/2 edges. The main abscissa for 2-connected graphs is equal to 
2.26 (see [3]), and for 3-connected graphs is equal to (7 + V / 7)/4 ~ 2.41 (see Theorem 16. 10|) . 




Figure 2: Cumulative degree distribution for 2-connected planar graphs with [in edges, [i £ (1, 3) 
and k = 2, . . . , 10 (left); and 3-connected planar graphs with \xn edges, fi € (3/2,3) and k = 
3,..., 10 (right). 



8 Concluding remarks 

Fusy [15] has designed a very efficient algorithm for generating random planar graphs uniformly 
at random. He has performed extensive experiments on the degree distribution on planar graphs 
with 10000 vertices and his experimental results, which he has very kindly shared with us, fit very 
well with the constants in Table [T] 
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Appendix 

The coefficients Ti(x,w) 

Let u stand for u(x, t(x)), which is the solution of 

_ (l + u)(3u-l) 3 
X ~ UU, ' 

The coefficients Ti{x,w) of the singular expansion of T(x,z,w) in Lemma 16.51 are given by 

(3u - lfw / P . N r- 

T ° {X > W) = 27648(3^ + 2«-l + «,)(« + IK " ( " + 1 " w)VP 

m / , (3m-1) 6 M; f P 2 ,0 P 2l l 

T 2 (x,w) - 



82944(3u 2 + 2u - 1 + w) 2 (u + l) 2 u 5 \ (9u + l) 2 y/p 
(3m- 1) 6 m;^2m(m+ l)(3u+ 1) ( fo ^ 2 2 P 3 



T 3 (x,w) = H '\ (3it- 1) 2 u>-9m 2 - lOu - 1 + 

V ; 373248(u+l) 3 u 6 V (u+l-w)\/P 

where 

P={u+1- w)(-(3u - l) 2 w + 81u 3 + 99u 2 + 19u + 1) 

and 

P = (27u 2 + 6u + l)w 2 + (-126u 3 - 150u 2 - 26m - 2)w + 81m 4 + 180m 3 + 118u 2 + 2Qu + 1 

P 2 = (1458u 5 + 3807u 4 + 900u 3 + 114u 2 - 6it - l)w 3 + 

(6561u 7 + 20898u 6 + 8532u 5 - 7281u 4 - 1635u 3 - 132u 2 + 30u + 3)w 2 + 

(-3645u 8 - 30942u 7 - 46494u 6 - 13230u 5 + 7536u 4 + 1590u 3 - 18u 2 - 42u - 3)w+ 

13122u 9 + 47385u 8 + 61560u 7 + 30708m 6 - 228u 5 - 4530w 4 - 872m 3 + 36u 2 + 18m + 1 

P2,i = (-54u 4 - 45u 3 + 57u 2 - 15u + l)w 4 + 

(-243u 6 + 27u 5 + 1278u 4 + 858u 3 - lllu 2 + 35u - 4)w 3 + 
(1944m 7 + 6507u 6 + 5553u 5 - 576u 4 - 1530m 3 + 15m 2 - 15u + 6)w 2 + 
(-1215u 8 - 6561u 7 - 11439u 6 - 7005u 5 + 231u 4 + 1229w 3 + 75u 2 - 15u - i)w 
+ 1458u 9 + 6561u 8 + 11376u 7 + 8988u 6 + 2388u 5 - 794u 4 - 512u 3 - 36m 2 + lOu + 1 
P 3 = -(3u - ifw 3 + (162m 4 + 135u 3 - 27m 2 - 3u - 3)w 2 

+ (8lu 5 + 243u 4 + 270m 3 + 138u 2 + 33m + 3)w + -81u 5 - 261m 4 - 298u 3 - 138u 2 - 21u - 1 

The coefficients Di(y,w) 

We proceed to give the coefficients Di(y, w) of the singular expansion of D(x, y, w) in Lemma 16.81 
Let t = t(y), for y e (0, oo), be the unique solution in (0, 1) of 

(l-2i) / t 2 (l-t)(l8 + 36t + 5t 2 ) 

y ~ (i + 3i)(i-*) oxp v W+W+^W+W 

Then, Dq(ij, w) is the solution of 



1 + D a = (1 + yw) exp 



( VS(D (t-l)+t) 
^4(3t+l)(D + l) 
D 2 (t 4 - 12t 2 + 20i - 9) + D (2t 4 + 6t 3 - 6t 2 + lOt - 12) + t 4 + 6t 3 + 9t 2 
4(t + 3)(A, + l)(3*+l) 
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where S is given by 

S = (D t -D Q + t)(D Q (t - l) 3 + *(* + 3) 2 ). 



The remaining coefficients D 2 (y,w), D 3 (y,w) are given in terms of Dq and *, 



(17* 5 + 237* 4 + 1155t 3 + 2527* 2 + 1808* + 400)(S 2 ,3 + S 2A VS) 



S-2,1 = -Dl{t - 1) 4 (* + 3)(lli 5 + 102* 4 + 411* 3 + 588* 2 + 352t + 72) 

- D a t(t - 1)(* + 3)(22* 7 + 231* 6 + 1059* 5 + 2277* 4 + 2995i 3 + 3272* 2 + 2000* + 432) 

- * 2 (* + 3) 3 (11* 5 + 85* 4 + 252i 3 + 108* 2 - 48* - 24) 

52,2 = D (t - 1)(11* 7 + 124* 6 + 582i 5 + 968* 4 - 977* 3 - 4828* 2 - 4112i - 984) 

+ *(* + 3) 2 (11< 5 + 85* 4 + 252i 3 + 108* 2 - 48* - 24) 
S^s = -(* + 3 )(A)i -D + t){D 2 {t - l) 4 + 2D {t - 1)(* 3 - * 2 + 5* - 1) + *(* 3 -'it- 14) 
S 2A = Dl(t 2 + 2t- 9)(< - l) 2 + D (2t 4 - 12t 2 + 80* - 6) + t(t 3 - 3* + 50) 

_ 24ft + 3)(£>p + l) 2 (t - l)t 2 (t + 1) 2 S 3 3 ^(S 3 , 2 - S 3 , 3 (D t -D + t)Vs) 
3 ~ 2 (D o t -D Q + t)(S 3A VS -(t + 3)(D * -D + t)S 3fi ) 

(3 = 3t(l + *)(17* 5 + 237* 4 + 1155* 3 + 2527* 2 + 1808* + 400) 

5 3 .1 = -5* 5 + 6* 4 + 135* 3 + 664* 2 + 592i + 144 

53.2 = Dl(Slt n + 135* 10 - 828i 9 - 180* 8 + 1982* 7 + 1090* 6 - 5196* 5 

+ 2108* 4 + 2425* 3 - 1617* 2 - 256* + 256) 
+ £>§(243i u + 1313i 10 + 1681i 9 - 51* 8 - 5269* 7 - 7325< 6 + 2571i 5 

+ 10271* 4 + 1846* 3 - 3888* 2 - 1392*) 
+ L>o(243i n + 2221* 10 + 8135i 9 + 15609* 8 + 12953* 7 - 3929* 6 - 12627* 5 

- 13293* 4 - 7680* 3 - 1632* 2 ) 
+ 81* 11 + 1043* 10 + 5626* 9 + 16806* 8 + 30165* 7 + 30663* 6 + 13344* 5 + 1008* 4 - 432* 3 
#3,3 = A)(81* 8 + 378* 7 + 63* 6 - 1044* 5 + 1087* 4 - 646* 3 - 687* 2 + 512* + 256) 
+ 81i 8 + 800* 7 + 3226* 6 + 7128* 5 + 8781* 4 + 4320* 3 + 384< 2 - 144* 

53.4 = £>o(* 4 - 12* 2 + 20* - 9) + A)(2* 4 - 12* 2 + 80* - 6) + * 4 - 3* 2 + 50* 

53.5 = D%(t 4 - 4* 3 + 6* 2 - 4* + 1) + D (2* 4 - 4* 3 + 12* 2 - 12* + 2) + * 4 - 3* 2 - 14* 

The coefficients Bi(y,w) 

Consider the singular expansion of the solution to the integral of Theorem 16. 31 



and define Iij (y, w) as the coefficients of the singular expansion of IiZ 1 in terms of X = yl — x/R(y), 
when replacing w by D(x,y,w)/E(x,y,w). That is, 




dt = I (x, w) + I 2 (x, w)Z 2 + h{x, w)Z 3 + 0{Z 4 ) 



Ii(x,w)Z 



ko(y^ w ) + ha{v^)x 2 + i h3 {y 7 w)x 3 



+ 0{X 4 ) 



These coefficients j are given by 
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/o '° = + 2048t6 ( 4(3< ~ 1){t + 1)3 l0gW + 8(3t ~ 1)(< + 1)3 l0g( * + 1} 

+ 8(3i 4 + 6i 2 -l)log(2) - 2(t- l) 3 (3t + l)log(A) - 2(3i 4 + 6t 2 - 1) log(B) 

+ {t- l)(D a (t 3 - 3t 2 + 3t - 1) + t 3 + At 2 + t + 2)VS 

- f —^(D 2 (t e - 2t 5 + t 4 - At 3 + lit 2 - lOt + 3) 
t + 3 y uv ; 

+ D a {2t 6 + 8t 5 - 10t 4 - 32t 3 + 46i 2 - 8t - 6) + t 6 + 10t 5 + 34i 4 + Ut 3 + 21t 2 + 18t 



I0.2 = (3t+ 5 ll^ (4(3* 4 - 4i 3 + 6t 2 1) log(2) + 2(3i 4 + 6t 2 1) log(t)+ 

4(3i 4 + 6t 2 - 1) log(< + 1) + (-3i 4 + 8t 3 - 6t 2 + 1) log(3A) + {-3t 4 - 8t 3 - 6t 2 + 1) log(B) 
D 2 {t - l) 5 (R Q ,aVS + R ,i) - ^^(R 0A VS + R , 5 ) - 



(t + 3) [Ro.6^S + R ,3^ 

j (3f + l)t 2 (t - l) 8 (i?o, Vff + i?p, 1 )( J 3 /? 5/2 + D a 3/2 R ,2) 

0,3 512(t + 3)i 8 /3 5 / 2 (i?o,6^ + i?o, 3 ) 

(f-l) 6 i? 2 , (i?2,lV£+i?2, 2 ) 

/2,2 ~~ 3072/?i 6 (t+ 1)(£) + 1) 

(i + 1) 2 (1 + 2t) 2 a 3 / 2 (3i + l)(t - l) 7 ((t + 3) 2 i? 2il ^ - fl 2 , 2 ) 
123 = 16t 5 f3 5 / 2 {l + D ) 

V3(t - l) 6 i?^o 2 (fls.oVS - R 3A (D t -D + t)) 
/3 ' 3 = 2304V3TTTt 5 /3 3 / 2 (t + l) 3 / 2 (^ t - Do + *) 
where the expressions A, B, and polynomials are given by 

A = D (5t 3 - 3t 2 - t - 1) + 5t 3 + 6t 2 + 5i + (3t + 1)VS 

B = D (t 3 - 3t 2 + 3t - 1) + t 3 + 2t 2 + 5t+(t- 1)VS 
Ro, = 3D 2 (t - l) 2 - D (7t - 3) - t(t + 3) 
Ro,i = 3Z?o(t - l) 4 - D' 2 (t - l)(3t 3 - t 2 + 25t - 3) 

+ D t{t 3 + 8t 2 + 21t - 14) + (t + 3) 2 t 2 
R ,2 = 128t(3t + l)(t - 1)(1 + 2t) 2 (t + 3) 2 (t + l) 2 

Ro,3 = D 2 (t - l) 4 + 2D (t - l)(t 3 + t 2 + 3t - 1) + t 4 + 4t 3 + 7t 2 + 2t + 2 

R 4 = 3D 3 (t - l) 5 (51t 8 + 1081t 7 + 8422t 6 + 31914t 5 + 59639t 4 + 42461t 3 + 7584t 2 - 2832i - 864) 

- D 2 (t - l) 3 (153t 10 + 3204t 9 + 29055* 8 + 146710i 7 + 43295U 6 

+ 717528i 5 + 561457t 4 + 208750t 3 + 47040< 2 + 13248t + 2592) 

- D (t + 3) 2 t(408t 10 + 6177t 9 + 34003i 8 + 92097i 7 + 122523i 6 

+ 126075t 5 + 145777t 4 + 82707t 3 - 1543< 2 - 15088t - 3312) 
3t(t - l)(t + 3) 2 (2t 4 + 3t 3 - 2t 2 + 3t + 2)(400 + 1808^ + 2527< 2 + 1155i 3 + 237t 4 + 17t 5 ) 
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i? 5 = 3D 4 (t - l) 7 (51t 8 + 108U 7 + 8422* 6 + 31914* 5 + 59639t 4 + 42461* 3 + 7584* 2 - 2832i - 864) 
+ 2D 3 1 (t - l) 4 (249i 10 + 3333< 9 + 22417* 8 + 105245* 7 + 339675* 6 + 
664087* 5 + 513315* 4 + 127943* 3 - 6936< 2 - 1152* + 1296) 

- Dlt(t - 1) 2 (357* 12 + 7089* 11 + 58637* 10 + 273500* 9 + 828314* 8 + 1886278* 7 + 3638786* 6 

+ 5441836* 5 + 4731 121t 4 + 1945329* 3 + 179665* 2 - 96240* - 20304) 

- 2D t(t + 3) 2 (51i 12 + 849* 11 + 6580* 10 + 33465* 9 + 115887* 8 + 253743* 7 + 285517* 6 

+ 148083* 5 + 130634* 4 + 141380* 3 + 59715* 2 + 4944* - 1200) 
+ 3* 2 (* - 1)(* + 3) 3 (2* 4 + 3* 3 - 2* 2 + 3* + 2) (400 + 1808* + 2527* 2 + 1155* 3 + 237i 4 + 17* 5 ) 
Ro,6 = D Q (t~l) 2 + t 2 + t + 2 

#2,0 = 3(3* + 1)0 + l)(-5* 5 + 6* 4 + 135i 3 + 664* 2 + 592* + 144) 
#2,1 = £» 2 (3t 3 - 12* 2 + 7* + 2) + D (6t 3 - 3t 2 + *) + 3t 3 + 9t 2 
#2,2 = D 3 (3t 7 - 47* 5 - 18< 4 + 21* 3 + 164* 2 - 105* - 18) 
+ D 2 (9t 7 + 36i 6 - 19* 5 - 168* 4 - 165* 3 + 292i 2 + 15*) 

+ A)(9* 7 + 72* 6 + 190* 5 + 156* 4 - 63* 3 - 108* 2 ) + 3* 7 + 36* 6 + 162* 5 + 324* 4 + 243* 3 
#3,0 = Dl{t 4 - 2t 3 + 2t - 1) + A)(2* 4 + 4t 3 - 2* 2 - 4*) + * 4 + 6* 3 + 9t 2 
#3,1 = Dl(t 5 - 3* 4 + 2* 3 + 2* 2 - 3* + 1) + D (2t 5 + 6< 4 - 2* 3 - 6* 2 ) + * 5 + 9* 4 + 27* 3 + 27* 2 

Finally, we write the coefficients Bi(y,w) in Lemma 16.91 in terms of the I+ j, Di, the radius of 
convergence R = R(y), and y and t. 

B = ^3 (-8 log(2)(3* 4 + 6< 2 - 1) - 8 log(* + 1)(3* - 1)(* + l) 3 - 41og(*)(3* - 1)(* + l) 3 
+ 21og(A)(* - 1) 3 (3* + 1) + 21og( J B)(3< 4 + 24* 3 + 6* 2 - 1) 
+ VS(t - 1)(A)(* 3 - 3t 2 + 3* - 1) + * 3 - 8* 2 + * - 2) 
Do 



t + 3 



(D (t - 1) 5 (* 2 + 2t - 9) + 2(* - 1) 3 )(* 4 + 60* + 3) + (f + 3) 2 (* - 1)(* 3 - 8* 2 + * - 2)*)) 

_ RDq(Dq(R 2 E 2 + RE 2 ) - 2(1 + REp)) 1 
B2 2{l + REo) 2 + # ( °'° + °' 2 + 2 ' 2 > 

+ (log(l + D ) - log(l + yw) - f^° ) 0- + D - D2)R 

R 2 D (2D 3 E%R + 2D 3 E + E 3 D Q ) 1 
B 3 = 2(k# + l) 2 + RD3{hg{1 + yw) ~ log(jDo + 1)} + ^ (/o ' 3 + /2 ' 3 + /3 ' 3) 
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